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A  NEW  TECHNIQUE  FOR  THE  DETERMINATION 
OF  HEATS  OF  SUBLIMATION  OF  SLIGHTLY  VOLATILE  SUBSTANCES 


By  O.  C.  Simpson  and  R.  J.  Thorn 


A  technique  is  described  which  enables  one  to  measure  relative  vapor  pressures  by  measuring 
the  rate  of  change  of  transmissivity  of  thin  deposits  as  they  are  formed  by  condensation  from  the 
vapor.  From  the  temperature  dependence  of  this  rate  the  heat  of  sublimation  of  slightly  volatile  sub¬ 
stances  can  be  determined.  In  the  case  of  graphite,  deposits  of  the  order  of  0.1  pg  per  cm2  can  be 
detected. 


INTRODUCTION 

Two  general  methods  of  determining  the  heat  of  vaporization  or  sublimation  are  discussed.  One 
method  is  based  upon  the  measurement  of  the  slope  of  the  plot  of  log-^p  against  the  reciprocal  abso¬ 
lute  temperature  and  requires  only  relative  vapor  pressure  measurements  as  a  function  of  tempera¬ 
ture;  the  other,  based  upon  the  calculation  of  the  absolute  entropy  difference  of  the  vapor  and  the 
condensed  phase,  requires  the  measurement  of  the  absolute  vapor  pressure  at  a  known  temperature 
plus  a  knowledge  of  the  pertinent  specific  heat  data.  Of  these  methods  the  second  is  usually  the  more 
accurate,  but  the  first  one  presents  certain  advantages  since  factors  such  as  the  vaporization  coef¬ 
ficient  must  be  rather  strongly  temperature  dependent  to  cause  large  error  and  since  only  relative 
vapor  pressure  measurements  are  required.  Also  no  specific  heat  data  are  needed. 

One  of  the  chief  difficulties  in  the  application  of  either  method  to  substances  having  low  vapor 
pressures,  graphite  for  example,  is  that  ordinary  techniques  for  determining  the  small  rates  of  evap¬ 
oration  require  long  times.  The  technique  described  in  this  report  enables  one  to  measure  extremely 
small  rates  of  evaporation  (actually  condensation)  and  hence  enables  one  to  employ  the  first  method 
without  employing  long  times.  The  same  technique  can  also  be  used  in  the  second  method  provided  a 
proper  calibration  is  made. 

The  authors  are  employing  this  technique  in  conjunction  with  the  first  method  to  determine  the 
heat  of  sublimation  of  graphite.  As  yet  the  investigation  has  not  been  completed  to  the  authors1  satis¬ 
faction.  However,  others  may  have  use  for  this  technique  and  since  it  will  be  some  time  before  the 
investigation  is  concluded,  the  experimental  method  is  described  here. 


THEORY 

The  technique  presented  herein  is  essentially  as  follows.  A  fraction  of  the  material  evaporating 
from  the  top  surface  of  a  cylinder,  heated  in  an  evacuated  space,  is  collected  on  a  circular  target  placed 
above  and  parallel  to  the  surface.  The  time  rate  of  deposition  of  the  material  is  determined  by  mea¬ 
suring  the  time  rate  of  change  of  the  transmissivity  of  the  target.  The  light  from  the  incandescent 
material  serves  as  a  light  source  for  this  measurement.  At  feasible  temperatures  for  substances  in 
which  we  are  interested,  graphite  for  example,  the  rate  of  deposition  is  so  small  that  ordinary  methods 
such  as  weighing  do  not  permit  its  determination  within  reasonable  times.  However,  the  use  of  an 
optical  pyrometer  constitutes  a  very  sensitive  method  for  the  determination  of  this  rate.  From  the 
dependency  of  the  deposition  rate  on  temperature  it  is  possible  to  calculate  relative  vapor  pressures 
and  therefrom  the  heat  of  sublimation.  The  experimental  arrangement  is  represented  schematically 
in  Figure  1. 
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Figure  1.  Experimental  arrangement  used  to  determine  rate  of  evaporation  of  solids. 

If  light  of  intensity  IQ  and  wavelength  X  emitted  from  the  top  surface  of  the  cylinder  at  some  high 
temperature  T°K  is  intercepted  by  the  layer  of  the  material  deposited  on  the  quartz  target,  then  the 
decreased  intensity  I  of  light  entering  the  optical  pyrometer  is  given  in  terms  of  the  thickness,  x,  of  the 
deposit  by  the  relation 

I  =  I0(l  -  rjHl  -  r2)(l  -  r3)  e_kx  ‘Vq 


In  this  equation  r^,  rg,  and  r^  are  the  reflectivities  of  the  deposit- vacuum,  the  deposit-quartz,  and  the 
quartz-vacuum  interfaces  respectively;  k  and  kq  are  the  absorption  coefficients  of  the  material  and 
quartz;  and  Xq  is  the  thickness  of  the  quartz.  More  correctly  there  may  be  additional  terms  due  to  the 
window  on  the  apparatus  and  a  prism.  These  terms,  however,  can  be  omitted  from  this  discussion, 
although  they  must  be  considered  when  it  is  necessary  to  measure  the  true  temperature  of  the  cylinder. 
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If  we  assume  that  the  spectral  intensity  I0  emitted  from  a  nonblack  body  with  a  spectral  emis- 
sivity  e^  is  given  by  Wien’s  law,*  we  can  write 


I0  =  exc1X'Vc2/XT 


(2a) 


Corresponding  to  the  decreased  intensity  I  there  is  some  apparent  absolute  temperature  Ta°K  such 


that 


I  =  e 


xc1x',5e_c2ATa 


(2b) 


If  equation  2b  is  divided  by  equation  2a  and  the  expression  for  I/I0  so  obtained  is  substituted  into 
equation  1,  then  an  equation  relating  the  apparent  temperature  Ta  to  the  thickness  and  reflectivities 
is  obtained.  That  is, 


T~a  =  T  +  ^  ^  +  ^  ^  ^  1(1  "  riK1  “  r2)(1  ~  r3>l  <3> 


If  the  thickness  x  of  the  material  deposited  on  the  quartz  is  increasing  with  time,  due  to  molecules 
condensing  on  the  target  (in  this  discussion  we  are  assuming  that  the  vapor  is  composed  of  a  single 
molecular  species),  the  dependency  of  Ta  on  time  can  be  found  from  the  elementary  kinetic  theory 
relations.  The  result  is 


1  1  t  Xkqxq  X 

Ta  T  C£  Cg 


In  [(1  -  ri)(l  -  rjjHl  -  r3)] 


XkFAs 
+  c2DAd 


|/H  J2_ 

y  27tr  yr 


t 


(4a) 


In  this  expression  F  is  the  fraction  of  the  molecules  evaporating  which  condenses  on  an  area  that 
intercepts  the  light  entering  the  pyrometer;  As  is  the  area  of  the  top  surface  of  the  cylinder,  A^  is  the 
area  of  the  deposit  through  which  light  passes  to  the  pyrometer;  D  is  the  density  of  the  deposit;  M  is 
its  molecular  weight;  R  is  the  molal  gas  constant;  p  is  the  vapor  pressure  of  the  material;  and  t  is  the 
time.  In  the  derivation  of  this  equation  there  is  the  inherent  fact  that  when  t  =  0,  r^  is  zero  and  r2  =  r3, 
since  no  deposit  is  present  on  the  target.  For  very  small  values  of  t,  r2  becomes  different  from  r3  and 
rj  becomes  greater  than  zero,  but  neither  is  necessarily  independent  of  t  (i.e.,  x), because  the  re¬ 
flection  of  light  takes  place  not  only  from  the  outer  atomic  layer  but  also  from  a  finite  number  of 
atomic  layers.  However,  for  sufficiently  large  values  of  t  these  reflectivities  become  constant  since 
the  deposit  is  then  sufficiently  thick  to  permit  the  maximum  reflection  characteristic  of  the  interfaces. 
Hence  for  values  of  t  greater  than  this  minimum  and  for  constant  values  of  the  temperature  of  the 
cylinder,  equation  4a  assumes  the  form 


i  *  Aipt 

=  A0  +  (4b) 

This  relation  shows  that  the  reciprocal  of  the  apparent  absolute  temperature  is,  under  the  conditions 
already  specified,  a  linear  function  of  time  with  a  slope, 

a  =  XkFAg  P  (5) 

C2DAd  '2irR 

Since  the  Clausius-Clapeyron  equation, 

dT  T4V  K  1 


*  For  X  =  0.665  and  T  =  2500°K,  eC^  =  5.5  x  10^.  Hence  one  is  justified  in  using  Wien’s  law 
rather  than  Planck’s.  Whether  Wien’s  law  is  strictly  applicable  to  nonblack  bodies,  however,  is  another 
question  which  will  not  be  discussed. 
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can  be  expressed  in  a  form  involving  the  derivative  of  the  logarithm  of  the  vapor  pressure,  determi¬ 
nation  of  the  temperature  rate  of  change  of  a  quantity  that  is  proportional  to  the  vapor  pressure  such 
as  the  slope  a  in  the  present  case  is  sufficient  for  the  determination  of  the  heat  of  sublimation. 
Equation  6  can  be  expressed  in  the  desired  form  by  neglecting  the  volume  of  the  solid  phase  com¬ 
pared  with  the  gaseous  phase  and  assuming  the  latter  to  be  a  perfect  gas.  If  the  value  of  p  given  by 
equation  5  is  then  substituted  in  this  relation,  we  find 


d  In  (a  VT)  __AH 


(7) 


Hence  the  heat  of  sublimation  is  equal  to  4.575  times  the  slope  of  the  plot  of  log1Q(a  j/T)  vs  1/T. 

Provided  the  absorption  coefficient  can  be  evaluated,  then  the  technique  described  herein  can  be 
used  to  determine  the  absolute  vapor  pressure  of  the  material  at  a  given  temperature  (see  equation  5). 
From  this  the  heat  of  sublimation  can  be  calculated  according  to  statistical  mechanical  relations.* 


♦ 


EXPERIMENTAL  EXAMINATION  OF  TECHNIQUE 


Typical  examples  of  the  plot  of  the  reciprocal  apparent  absolute  temperature  Ta  measured  through 
graphite  deposits  as  a  function  of  the  time  are  given  in  Figure  2.  These  curves  are  characterized  by  a 
slight  negative  curvature  in  the  initial  region  which  becomes  less  with  increasing  time  until  the  curve 
finally  becomes  linear.  The  length  of  the  nonlinear  region  depends  upon  the  temperature  of  the  cylinder, 
increasing  with  decreasing  temperature.  This  initial  curvature  is  to  be  expected  according  to  equation 
4a,  since  during  this  period  the  reflectivities  rj  and  T2  are  time  dependent  and  are  changing  in  such  a 
manner  as  to  make 

d  In  [(1  -  ti)(1  -  r2)] 
dt 


a  negative  quantity  that  decreases  numerically  to  zero  at  r^  and  r2  constant.  Hence  1/Ta  vs  time  has  a 
negative  curvature  until  r1  and  r2  become  constant.  The  possibility  of  this  nonlinear  region  having  an 
effect  on  the  calculated  slope  a  can  be  excluded  in  two  ways.  In  one  procedure  the  slope  of  only  the 
linear  region  is  calculated,  this  region  being  defined  as  follows.  If  tj  is  the  minimum  time,  which  for 
all  practical  purposes  is  necessary  for  the  attainment  of  the  linear  region  at  some  temperature  Tj, 
then  the  equivalent  time  t2  necessary  at  some  other  temperature  T2  can  be  found  to  be 

(8, 


Since  the  linear  region  is  readily  discernible  at  the  high  temperature  end  of  the  permissible  range,  one 
of  these  runs  can  be  selected  arbitrarily  for  estimating  tj.  The  corresponding  t’s  for  other  tempera¬ 
tures  can  then  be  calculated  with  the  aid  of  this  equation  provided  a  rough  estimate  of  A  H  has  been 
made.  Since  accurate  values  of  the  t’s  are  not  required,  AH  can  be  estimated  by  considering  1/Ta  vs  t 
to  be  linear  for  all  values  of  t.  In  essence  this  procedure  is  equivalent  to  selecting  regions  of  equiva¬ 
lent  thickness.  Another  procedure  of  eliminating  the  initial  region  is  by  an  initial  deposit  at  a  rela¬ 
tively  high  temperature  before  measurements  of  Ta  vs  time  are  made.  When  this  procedure  is  used 
1/Ta  is  linear  in  time. 


The  deposits  whose  transmissivities  are  represented  in  Figure  2  range  in  total  weight  from 
23.7  (ig  at  2252°K  to  143  pg  at  2433°K  for  an  area  of  about  3  cm2.  Since  one  can  readily  read  a  tem¬ 
perature  difference  of  3  °  with  an  optical  pyrometer,  the  limit  of  detection  of  the  presence  of  the 
deposit  is  0.11  /ig  per  cm2.  This  corresponds  to  a  thickness  of  about  5  x  10“ 8  or  about  three  or  four 


T 


* K.  K.  Kelley,  U.  S.  Bur.  Mines,  Bull.  383  (1935). 
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atomic  layers.  This  conclusion  agrees  with  that  of  Estermann  and  Stern,*  who  found  that  deposits  of 
two  or  three  molecules  in  thickness  are  visible  because  of  aggregation.  The  thickness  of  the  graphite 
at  which  the  curvature  becomes  zero  or  at  which  the  reflectivities  become  constant  is  about  3.6  x  10"® 
cm  or  about  250  atomic  layers. 


Figure  2.  Relative  rates  of  deposition  of  graphite. 


♦Estermann  and  O.  Stern,  Z.  Physikal  Chem.  106:399  (1923). 
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POSSIBLE  ERRORS 

The  possible  sources  of  error,  which  may  be  considered  as  being  inherent  in  the  transmissivity 
technique,  are  (1)  the  dependency  of  the  absorption  coefficient  on  thickness  and  (2)  on  temperature, 

(3)  the  deviation  of  nonblack  body  radiation  from  Wien's  law,  (4)  thermal  expansion  of  the  deposit,  and 
(5)  the  radiation  from  the  deposit. 

In  general  the  effect  of  factors  (1)  and  (3)  if  operative  would  be  to  make  1/Ta  other  than  a  linear 
function  of  time.  Since  we  have  observed  a  linear  dependence  with  graphite  aside  from  the  initial 
effect  of  reflectivities,  we  conclude  that  these  factors  cause  no  error  in  this  case  unless  a  fortuitous 
combination  of  them  causes  effective  cancellation.  It  is  doubtful  that  the  absorption  coefficient  of 
graphite  changes  sufficiently  with  temperature  to  cause  any  error  since  the  transmissivity  curve  of 
thin  graphite  deposits  as  a  function  of  wavelength  shows  no  absorption  peaks  in  the  neighborhood  of 
0.665  p.  The  thermal  err^ansion  of  graphite  would  cause  log  a  j/T  to  change  by  less  than  0.02%  for  the 
temperature  range  given  in  Figure  2.  If  the  deposit  on  the  quartz  target  emits  radiant  energy  and  if 
this  energy  were  of  the  same  order  as  that  from  the  cylinder  transmitted  through  the  deposit,  then 

values  of  a  will  be  in  error.  However,  the  ratio  of  these  intensities  is  e  ■c2ATd/re-c2AT 

in  which 

T(j  is  the  temperature  of  the  deposit  and  t  its  transmissivity.  Crude  measurements  showed  that 
r  =  0.01  (for  the  thickest  deposit)  and  T^^  1000 °C.  Hence  for  the  ratio  we  obtain  8  x  10“^  at  T  =  2200°K 
showing  that  the  amount  of  radiant  energy  emitted  by  the  deposit  is  negligible  in  our  experiments  com¬ 
pared  with  that  transmitted  from  the  cylinder. 

The  sources  of  error,  which  are  not  necessarily  inherent  in  the  transmissivity  method  but  which 
may  be  present  under  the  conditions  of  the  experimental  procedure,  are  (1)  temperature  gradients  in 
the  material,  (2)  variation  of  the  vaporization  coefficient  with  temperature,  and  (3)  variation  of  the 
condensation  coefficient  with  temperature.  For  metals  temperature  gradients  will  usually  cause  no 
trouble  since  their  thermal  conductivities  are  high  and  temperatures  above  2000 °K  are  not  usually 
required.  However,  at  temperatures  above  2000  °C  these  gradients  become  appreciable,  especially  in 
the  case  of  graphite,  and  they  must  be  taken  into  consideration  by  calculation  of  an  “average”  tempera¬ 
ture.  One  of  the  best  ways  of  eliminating  the  variation  of  the  vaporization  and  condensation  coefficients 
with  temperature  is  to  employ  conditions  such  that  they  are  unity.  With  the  vaporization  coefficient  this 
can  be  assured  if  the  material  is  vaporized  in  a  small  furnace  cavity  and  permitted  to  effuse  through 
an  orifice  in  the  furnace.  In  all  probability  the  vaporization  coefficient  will  also  be  unity  if  the  surface 
of  the  material  is  chemically  clean.  With  the  condensation  coefficient  a  value  of  unity  can  be  approached 
by  employing  a  water-cooled  target. 
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